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A two-photon transition is used to convert an arbitrary 
fraction of the 87 Rb atoms in a \F — 1, rrif = —1} condensate 
to the \F = 2, to/ = 1) state. Transferring the entire pop- 
ulation imposes a discontinuous change on the condensate's 
mean-field repulsion, which leaves a residual ringing in the 
condensate width. A calculation based on Gross-Pitaevskii 
theory agrees well with the observed behavior, and from the 
comparison we obtain the ratio of the intraspecies scattering 
lengths for the two states, aii,-i> /ai2,n=h062(12). 



The effects of interactions on dilute atomic Bose- 
Einstein condensates (BEC) jjj have been studied in sev- 
eral contexts, including density and momentum distribu- 
tions collective excitations ||], specific heatGj, speed 
of sound ||, and limited condensate number |(|7 There 
has been excellent quantitative agreement between these 
experiments and theory, made possible by both experi- 
mental advances and the fact that the interactions can 
be modeled relatively simply. Despite complicated in- 
teratomic potentials, the mean-field interaction in con- 
densates is well-characterized by a single parameter, the 
s-wave scattering length a. Previous quantitative exper- 
iments on interactions have all been done using single- 
component condensates with a constant value of a, but 
several authors have proposed using external optical or 
magnetic fields to shift the mean-field interaction by per- 
turbing the interatomic potential 

We present a method for the creation of condensate 
mixtures using radiofrequency (rf) and microwave fields. 
We are able to transfer abruptly a trapped condensate 
of one hyperfine state into a coherent superposition of 
two trapped hyperfine states, and then watch the subse- 
quent dynamical behavior. This approach makes possible 
a variety of two-species BEC studies. In this paper, we 
examine quantitatively one special case, in which all con- 
densate atoms are converted from one state to another. 
Since these two states have slightly different values of a, 
the sudden change in self-interaction gives rise to oscil- 
latory spatial behavior of the condensate wave function 
. The scattering length ratio can be extracted from a 
model using analytical equations of motion for the con- 



densate widths |T(^-|T2| based on Gross-Pitaevskii theory. 

The first demonstration of a binary mixture of conden- 
sates by Myatt et al. Q produced overlapping conden- 
sates of the 5S±/2 \F = 1, mf = —1) and \F = 2, m/ = 2) 
states of 87 Rb. The ratio of the magnetic moments of 
these states is 1:2, so the condensates experience dif- 
ferent potentials in a magnetic trap and are displaced 
unequally from the trap center by gravity. Due to an ac- 
cidental coincidence between the singlet and triplet scat- 
tering lengths of 87 Rb, collisional loss is reduced and any 
mixture of spin states will be relatively long-lived [[lj]-[l5) . 
Here, we use mixtures of |1, —1) and |2, 1) states, which 
possess several advantages. First, these two states have 
essentially identical magnetic moments, and hence feel 
identical confining potentials. Second, one can conve- 
niently and quickly change atoms from the |1, — 1) state 
to the 1 2, 1) state by a two-photon transition (microwave 
plus rf). Finally, we can selectively image the different 
components using appropriately tuned lasers. 

The apparatus uses a multiply loaded, double 
magneto-optical trap (MOT) scheme |ll|, consisting of 
two MOT regions connected by a 60 cm long, 1 cm di- 
ameter transfer tube. In the first MOT, 87 Rb atoms are 
collected from a background Rb vapor, produced by heat- 
ing a Rb getter source |llj . The trapping beams are then 
turned off and a near-resonant laser beam is pulsed on to 
push the atoms through the transfer tube into the second 
MOT. To prevent collisions with the walls, the outside 
of the transfer tube is lined with strips of permanent 
magnetic material which form a hexapole guiding field. 
The atoms exit the tube and are captured by the second 
MOT, which has a 70 second lifetime due to differential 
pumping. This process is repeated twice a second for 
~10 s until approximately 10 9 atoms have accumulated 
in the second MOT. 

The atoms are further cooled in an optical molasses 
and then optically pumped into the |1, — 1) state. These 
atoms are then transferred into a time-averaged, orbit- 
ing potential (TOP) magnetic trap @. The TOP trap 
consists of a quadrupole magnetic field with an axial gra- 
dient of ^250 G/cm and a uniform magnetic field which 
rotates at 1.8 kHz. The fraction of atoms captured from 
molasses in the TOP trap depends approximately linearly 
on the magnitude of the rotating field. We use the maxi- 
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mum rotating field achievable in our apparatus (49 G) to 
capture up to 50% of the atoms from the MOT, substan- 
tially enhancing the initial number as compared to our 
previous setup Q. As in (!]], the confined cloud is then 
magnetically compressed and evaporatively cooled to a 
temperature where there is no noticeable non-condensed 
fraction remaining, corresponding to a temperature less 
than 0.4 of the BEC transition temperature. Evaporative 
cooling requires about 30 seconds, and typically produces 
a condensate of 5 x 10 5 atoms. 

The magnetic moments of the |1, —1) and |2, 1) states 
are nearly the same to first order in magnetic field. 
Second-order dependence, the nuclear magnetic moment, 
and small effects due to the time varying nature of the 
TOP trap must be taken into account when calculating 
the exact trap potential Jl^,^]. The quadrupole and 
rotating magnetic fields are adiabatically changed after 
evaporation to make the trap potentials the same (to 
within 0.3%) for the two states. Evaporation takes place 
in a trap with radial frequency v r = 35 Hz; subsequent 
work occurs at v r = 17 Hz with a rotating field magni- 
tude of 5.7 G [y r = v z /VS for the TOP trap). 

The two-photon transition used to change the hyper- 
fine state is shown schematically in Fig. |]a. We apply 
a pulse of microwave radiation at a frequency slightly 
less than the ground state hyperfine splitting of 87 Rb 
(~6.8 GHz) along with a -2 MHz rf magnetic field. This 
connects the |1, —1) state to the |2, 1) state via an inter- 
mediate virtual state with a detuning of 2.2 MHz from the 
1 2,0) state . The two-photon Rabi frequency is 600 Hz, 
which is much faster than the characteristic frequency 
2v z = 96 Hz for the condensate to change shape. 

We are able to view either component of a mixture 
of |1, — 1) and |2,1) condensates or the combined den- 
sity distribution. Observation of the atoms is a destruc- 
tive measurement in which the TOP trap is suddenly 
turned off and the condensate is left to expand ballisti- 
cally for 22 ms. To image the combined density distri- 
bution of both states, a short pulse of "repump" light 
pumps the atoms from the |1, — 1) state into the F = 2 
manifold. About 100/is later, a a + circularly polarized 
probe beam 17 MHz detuned from the 5Si/ 2 F — 2 to 
5P 3 / 2 F' = 3 cycling transition is scattered by the atoms 
and the shadow is imaged onto a charge- coupled device 
camera. Imaging only the |2, 1) condensate uses the same 
procedure, except that the pulse of repump light is omit- 
ted. Atoms in the 1 1 , — 1 ) state are far (6.8 GHz) from 
resonance and invisible to the probe beam. Imaging only 
the |1,— 1) atoms is similar to viewing both species si- 
multaneously, but the |2, 1) atoms are first "blown away" 
from the imaging region by a 2 ms, 60 /iW/cm 2 pulse of 
F = 2 to F' = 3 light applied near the end of the ballistic 
expansion. This light has no effect on the |1, — 1) atoms, 
to which a subsequent pulse of repump light is applied, 
pumping them into the F = 2 manifold for probing. In 
Fig. [Tjb we show the population of each state as a func- 
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FIG. 1. (a). A diagram of the ground state hyperfine levels 
(F = 1,2) of 87 Rb shown with Zeeman splitting due to the 
presence of a magnetic field. The two-photon transition is 
driven between the |1,— 1) and |2, 1} states, (b) The Rabi 
oscillation of population between the 1, — 1) (open circles) 
and |2, 1) states (solid circles) as a function of the two-photon 
drive duration. The lines are fit to the data and show the 
expected sinusoidal oscillation. 



tion of the two-photon pulse length. The Rabi oscilla- 
tions have nearly 100% contrast, indicating that we are 
able to put any desired fraction into the |2, 1) state. 

We have investigated the dynamical behavior of the 
condensate after more than 99.5% of the atoms are trans- 
ferred from the |1, —1) state to the |2, 1) state. After the 
transfer, the condensate is allowed to evolve in the trap 
for a time interval T before the cloud is released from 
the trap and probed. The experiment is repeated for 
different values of T, and the axial and radial widths 
of the condensate are measured using a fit to the two- 
dimensional condensate image Fig. || illustrates the 
time dependence of both the axial and radial size of the 
atom cloud. Qualitatively, the data for both dimensions 
are consistent with a "compression" oscillation. The con- 
densate shrinks for very early times, indicating a weaker 
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FIG. 2. Oscillation in the width of the cloud in both the 
axial and radial direction due to the instantaneous change in 
scattering length. The widths are for condensates as a func- 
tion of free evolution time in units of the radial trap period 
(cdr 1 = 9.4 ms) , followed by 22 ms of ballistic expansion. Each 
point is the average of approximately 10 measurements. Note 
that the fractional change in width is quite small. The solid 
line is a fit of the model to the data, with only the amplitude 
of the oscillation and the initial size as free parameters. 

mean- field repulsion in the |2, 1) state. 

In order to describe the response to the change in scat- 
tering length, we use the Gross-Pitaevskii (GP) equation 
for the condensate wave function, with a time-dependent 
interaction term: 
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where V(r) is the confining potential, m is the Rb mass, 
and U(t) = 4Trh 2 a(t) / m. This has the form of a nonlinear 
Schrodinger equation, with a mean-field characterized by 
the s-wave scattering length a(t). For t < 0, a(t) = a%, 
the scattering length for |1, — 1) atoms on |1, — 1) atoms. 
For t > 0, a(t) = a 2 , the |2, 1) on \2, 1) scattering length. 
The condensate is formed in its ground-state in which 
the density profile n(r) = |<!>(r)| 2 is constant in time, 
and with the axial and radial widths (w z , w r ) determined 
by V(r), a\, and the number of atoms N. After a sud- 
den change in a(t), <&(r, f) is projected onto a coherent 
superposition of its new ground-state and collective vi- 
brational modes. Instead of solving the complete GP 
equation, one can use the Thomas-Fermi (TF) approx- 
imation, which corresponds to neglecting the quantum 
pressure term in the kinetic energy of the condensate. In 
this case, one can replace the GP equation with a pair 
of scalar equations of motion for the condensate widths 

,„ N2 fl5h 2 Na(t)\ l 
w t + (2^)V - ( „ B/2 „ 2 U ) = . (2) 



Here the widths w r = \J (r 2 ) and w z = \J (z 2 ) are the 
radial and axial rms radii, and the subscript i is to be 



replaced with either r or z for the respective widths. The 
TF approximation holds when N is large [p2| and, in this 
limit, the GP theory coincides with the hydrodynamic 
theory of superfluids Q . 

Time-dependent behavior is predicted by numerically 
integrating Eqs. (2). The initial conditions (t < 0) are 
that the widths are at the values determined by the 
ground-state solution, w r = w z = 0, and alt) = a%. 
After t = 0, a(t) — a 2 and the integration proceeds for 
time T. We model the removal of the trap by setting 
v r = v z = and continue the numerical integration for 
the 22 ms free expansion. The axial and radial widths 
subsequently expand with different speeds. In the TF ap- 
proximation, the ratio between the two expanded widths, 
averaged over T, can be shown to depend only on u r , v z 
and the expansion time [ [HI . In our case, the theory 
predicts w z /w r = 1.29, in good agreement with the mea- 
sured value w z /w r — 1.31. The oscillations of the widths 
correspond to a superposition of two m = modes. Since 
no angular momentum is imparted to the condensate by 
the change of scattering length, these are the only modes 
excited. The calculated mode frequencies are vjv r = 1.80 
and vjv r = 4.99, and turn out to be independent of the 
amplitude in the range considered here. The two modes 
contribute to the oscillations in z and r with a different 
phase and amplitude, the axial motion behaving mainly 
as the fast mode and the radial motion as the slow one. 

In each direction, the ratio of the amplitudes of the two 
modes is predicted to be constant over the range of possi- 
ble scattering lengths relevant to the experiment. Thus, 
only the initial size and an overall amplitude are used as 
adjustable parameters in comparison of theory with ex- 
periment. The solid lines in Fig. || show the theoretical 
prediction using the best fit value of the oscillation ampli- 
tude, which is related to the scattering length ratio ai/a 2 . 
The predicted oscillations agree remarkably well with the 
shape of the data, for both the frequency and the phase 
of the two modes. From the fit amplitude we obtain the 
ratio of the scattering lengths ai/a 2 = 1.062 ± 12, which 
is consistent with the ratio of 1.059lg obtained in a the- 
oretical calculation of binary collision parameters J|4|] . 

Analysis of the data is complicated slightly by an ob- 
served systematic dependence of the oscillation ampli- 
tude on the magnitude of the rf drive during the two- 
photon pulse. We believe that this effect is partly due to 
a small dressing of the atoms to the |1,0) or |2,0) state 
by the rf, and partly due to a coupling of the rf onto 
electronics controlling the trap potential. Both result in 
a change in strength of the confining potential during the 
two-photon pulse which returns to normal when the pulse 
is complete. After this impulsive perturbation, the BEC 
can be thought of as freely evolving with an initial "veloc- 
ity" in the width (w r ,w z ^ 0). This is in contrast to the 
discrete change in scattering length, in which the BEC 
width has an initial offset from the equilibrium width, 
but no initial velocity. The systematic is manifested as 
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FIG. 3. Dependence of the oscillation amplitude on the rf 
strength in the two-photon transition. The vertical axis rep- 
resents the fitted value of the fractional excitation amplitude, 
with error bars from uncertainty in the fit. The solid line 
shows extrapolation to zero rf amplitude. 



an initial offset in the phase of the resulting oscillation, 
which is indeed observed for large rf amplitudes. Our 
result for the ratio was obtained by extrapolating to zero 
rf amplitude with a quadratic fit as shown in Fig. ^. 

We have explored the dynamical response of a Bose- 
Einstein condensate to a sudden change in the interaction 
strength. The agreement with the model is excellent, and 
presents a new and precise way of measuring a ratio of 
scattering lengths. This also allows a test of the theory 
for pressure shifts in atomic clocks psj . Our method of 
creating two fully interpenetrating condensates will allow 
coherent control of the relative population and detailed 
examination of spatial and phase dynamics in the future. 
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